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ABSTRACT 

Herschel PACS and SPIRE images have been obtained over a 30'x30' area around the well-known carbon star CW Leo (IRC +10 216). An 
extended structure is found in an incomplete arc of ~22' diameter, which is cospatial with the termination shock due to interaction with the 
interstellar medium (ISM) as defined by Sahai & Chronopoulos from ultraviolet GALEX images. Fluxes are derived in the 70, 160, 250, 350, and 
550 /im bands in the region where the interaction with the ISM takes place, and this can be fitted with a modified black body with a temperature 
of 25±3 K. Using the published proper motion and radial velocity for the star, we derive a heliocentric space motion of 25.1 km s"'. Using 
the PACS and SPIRE data and the analytical formula of the bow shock structure, we infer a de-projected standoff distance of the bow shock of 
Rq = (8.0 ± 0.3) X 10" cm. We also derive a relative velocity of the star with respect to the ISM of (106.6 ± 8.7)/ y^isM km s"', where Uism is the 
number density of the local ISM. 
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1. Introduction 

Ever since the discovery paper by Becklin et al. (1969), the ob- 
ject IRC -HlO 216 (= AFGL 1381 = CW Leo) has spuiTed much 
interest. We now know that it is a carbon star in an advanced 
stage of stellar evolution on the asymptotic giant branch (AGB), 
pulsating and surrounded by an optically thick dust shell and 
large molecular circumstellar envelope (CSE). One aspect of 
study has been to constrain the properties of the CSE by answer- 
ing questions such as what is the mass-loss rate and how has it 
changed with time, what kind of chemistry takes place, and what 
is the geometry and structure of the CSE? 

The deep optical images taken by Mauron & Huggins (1999, 
2000), Mauron et al. (2003), and Leao et al. (2006) show that the 
dusty envelope is not smooth but consists of a series of arcs or 
incomplete shells. The average angular separation between the 
dust arcs suggests a timescale for the change in mass-loss rate 
of the order of 200-800 yr. The lack of kinematic information 
on the dust arcs precludes any firm conclusion about the true 
three-dimensional structure of the arcs or shells. 

From large-scale mapping at a relatively low angular resolu- 
tion of the CO J=l-0 emission, Fong et al. (2003) discovered a 
series of large molecular arcs or shells at radii of ~100" in the 
outer envelope. They attribute these multiple shells as "being the 
reverberations of a single Thermal Pulse erupting over 6000 yr 
ago." The timescale inferred from the spacing between these arcs 
is about 200-1000 yr. 

In the present paper, we discuss our new results on the outer 
shell of CW Leo from observations with the Herschel Space 
Observatory (Pilbratt et al. 2010) and their connection to the 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



results from GALaxy Evolution Explorer Space Observatory 
(GALEX) by Sahai & Chronopoulos (2010, hereafter SC). For 
the present analysis, we adopted a distance of d = 135pc and a 
mass-loss rate M- 2.2-10"'' Mq yr"' (Groenewegen et al. 1998), 
a gas expansion velocity of Vexp = 15.4 km s"', a radial ve- 
locity Vlsr = -25.5 km s"' (Groenewegen et al. 2002), cor- 
responding to Vheiio = -18.6km s"', and a proper motion (pm) 
Pa COS 6 = +26±6, ps = +4±6 mas/yr (Menten et al. 2006). 



2. Observations and data reduction 

The observations were carried out using the Photo AiTay Camera 
& Spectrometer (PACS, Poglitsch et al. 2010) and the Spectral 
and Photometric Imaging REceiver (SPIRE, Griffin et al. 2010), 
and are part of the Science Demonstration Phase observations 
of the Mass-loss of Evolved StarS (MESS) guaranteed time key 
program (Groenewegen et. al. 2010, in prep.), which is inves- 
tigating the dust and gas chemistry and the properties of CSEs 
around a large sample of post-main-sequence objects. 

For both instruments, the scan-map observing mode was 
used for a total sky coverage of 30'x30'. In this mode, the tele- 
scope is slewed at constant speed (20"/s for the PACS data, 30"/s 
for the SPIRE data) along parallel lines to cover the required area 
of the sky. For the PACS observations, two scan maps were taken 
with a scanning angle of 90° between the two to achieve the most 
homogeneous coverage. The observation identification numbers 
(obs.ID) for the two scans are: 1342186298 and 1342186299, 
with a total integration time of 8.35 hours. For the SPIRE ob- 
servations, a single scan was taken (obs.ID 1342186293) with a 
total integration time of 1.35 hours. All observations were taken 
on 25 Oct. 2009. 

The data reduction was performed using the Herschel 
Interactive Processing Environment (HIPE) version 2.3. We fol- 
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Fig. 1. Surface brightness maps in Jy/pixel. From left to right : PACS 160 jum image. SPIRE 250 yum image for which we overplotted 
the annuli segments between which we integrated the flux for the extended emission (inner annulus) and the sky (outer annulus). 
The white cross represents the center of the ellipse. SPIRE 250 jum image with overplotted in white the contour from the FUV map 
at 4.4x 10"^ mJyarcsec"^ limit. SPIRE 350 fim image. The FOV for all images is 23'x27'. Background sources were removed from 
all maps. North is up and east is to the left. 



lowed the basic data reduction steps as described in the PACS 
data reduction guide vl.2 and the SPIRE data users manual vl.O 
with the exception of the median filtering of the background. In 
the SPIRE pipeline, the entire scanline is used to calculate the 
median value. In the PACS pipeline, a running median filter is 
used with a certain "high pass filter width". A value of 450 was 
adopted that leads to a spatial scale of filtering approximately 
equal to that used for the SPIRE data, and is of the order of 30'. 
For both PACS and SPIRE, the central star was masked before 
applying the median filter The mask was a circle around the cen- 
tral star that approximated a contour at 3 cr of the background 
noise level. 

Near-UV (NUV) and fai--UV (FUV) GALEX data (Martin 
et al. 2005) were added to the analysis. The data were taken in 
2008 and ai-e centred around 1528A (FUV) and 2271A (NUV). 
We used the pipeline product retrieved from the GALEX archive. 

3. Analysis 

In the present paper, we do not comment on any structure close 
to the star previously detected in the optical and in CO. No sim- 
ilar structure is clearly visible in our data but a detailed analysis 
of the central structure requires a very detailed understanding 
and accurate subtraction of, or deconvolution with, the compli- 
cated PSF. This is beyond the scope of this paper 

The new result presented here is the extended emission in the 
form of an arc clearly seen at 160 yum, 250 /vm, and 350 /vm with 
a spatial scale as large as 22' (see Fig. [U. The arc is eastward 
of the central star with an almost circular curve. The central star 
is slightly to the east of the centre of the structure and the arc is 
slightly flattened in the easterly direction. It is interesting to note 
that this extended emission matches the position and the shape of 
the FUV extended emission (see Fig.[T]i. In the Herschel images, 
we do not see the patches of nebulosity seen in the western part 
of the FUV image (see Fig. 1 in SC). 

The ID intensity profile for each map (Fig. |2]i was con- 
structed by subtracting the total flux integrated within two suc- 
cessive apertures. To match the shape of the arc, elliptical aper- 
ture photometry was used. We only considered the area of the 
ellipse between 14° position angle (PA) and 167° PA (see Fig.[lJ 
to match the spatial scale of the extended emission. The ellipses 
are centred on o- = 146.955° and 6 = -h13.28° with a PA of -12° 
and a major-to-minor axis ratio taken to be 1.1. For the emis- 
sion from the sky, the inner minor axis is 13.3' and the outer 



minor axis is 13.8' (see Fig[T]i. The PACS fluxes were corrected 
following the numbers in the PACS Scan Map release note yield- 
ing an absolute flux calibration uncertainty of 10-20%. A sim- 
ilar correction was applied to the SPIRE fluxes following the 
correction factors given by Griffin et al. (2010) and Swinyard 
et al. (2010). The SPIRE absolute flux calibration uncertainty 
is ~15%. The dust shell is detected at 160 fim, 250 fj.m, and 
350 jum around 9.6' from the star for a maximal surface bright- 
ness of 0.067 mJy-arcsec"^ at 160 fim. There is an offset of ~20" 
between the FUV intensity peak and the far-IR intensity peaks, 
which suggests that they have different origins. 

Based on the location of the extended emission at 160 //m, 
the total flux in the arc was calculated within the segment of an 
elliptical annulus of 150" width, and inner minor axis of 8.5' be- 
tween PA 14° and 167° (see FiglTJ. The derived fluxes are listed 
in Table [T] The errors include the absolute calibration error and 
the error in the background estimation. We fitted the photometry 
of the extended emission with a function of the form By ■ A'^, ex- 
pected for a grain emissivity as Qabs ~ 'i^- The NUV and FUV 
data are plotted for reference, but these do not fit because the 
emission has a diff'erent physical origin (see below). The best fit 
solution infers a dust temperature of 25+3 K andjS = 1.6+0.4 
(see Fig. [3]). The value for /3 is in-between that expected for 
amorphous carbon (Rouleau & Martin 1991, y6= 1.1) and astro- 
nomical silicates (Volk & Kwok 1988,/?= 2.0), suggesting that 
the material in the bow shock is a mixture of C-rich material, 
lost by the star, and swept-up ISM material. 



4. Discussion 

4.1. Bow shock, thermal pulse, or both ? 

The most widely accepted explanation of large detached shells 
is mass-loss variation (e.g., Olofsson et al. 1990, Zijlstra et al. 
1992). AGB stars experience thermal pulses (TPs) during which 
intense mass loss ejections occur. A star can undergo several 
TPs, which would lead to the formation of concentric spherical 
detached shells (see Kerschbaum et al. 2010). 

Another explanation of detached envelopes is the interaction 
between the AGB wind and the interstellar medium (ISM) (e.g. 
Young et al. 1993, Martin et al. 2007, Ueta et al. 2006). In this 
scenario, the AGB wind is slowed down as it sweeps up material 
from the ISM. The piled up material forms a density enhance- 
ment that continue to expand due to the internal pressure. Shocks 
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Fig. 2. Intensity profiles as a function of the offset along the mi- 
nor axis of the ellipse increasing from west to east. The profile 
for each wavelength was normalised to the intensity at 558" and 
shifted up for clarity. From bottom to top with normalisation 
factor in mJy/arcsec^ and the vertical shift: FUV (0.065,0.0), 
160 jum (0.063,0.8), 250 ^im (0.023,1.1), 350 yum (0.008,1.4), 
500 /zm (0.004,1.8). The vertical dotted line indicates the posi- 
tion of the dust emission peak from the center of the ellipse. 



wether the dust shell we see is from an earlier TP that is now 
interacting with the ISM. 

Guelin et al. (1995) constrained the initial mass of CW Leo 
to be 3 < M < 5 Mq based on the isotopic ratios of 24,25,26]y|-g 
From the initial-final mass relation from Salaris et al. (2009), 
this implies a likely final mass (and essentially the current core 
mass) of 0.7-0.9 M©. From the core mass interpulse period rela- 
tion of Wagenhuber & Groenewegen (1998) for solar metallicity 
stars, this implies an interpulse period of 6 000-33 000 years. 

Scaling the mass loss and the wind velocity in SC with 
our assumed values, the flow timescale in the unshocked and 
shocked winds is 19 900 and 56 000 yr (very long because of 
the very low velocity of 1.2 km s"' in this region, see SC for 
details), for a total lower limit to the duration of the mass-loss 
phase of about 75 000 years. If TPs were to modulate the mass 
loss, as hypothesised for the origin for the detached shells dis- 
cussed earlier, one could expect at least one other TP to have 
occurred during the time it took the envelope to expand this far 
No obvious density enhancement is evident in the unshocked 
wind in the PACS and SPIRE images, which suggests that the 
interpulse period is at least 19 000 years. However, the dynam- 
ical and interpulse timescales are compatible, so it is possuble 
that the bow shock and dust emission are not only the result of 
a steady outflow interacting with the ISM, but might include the 
effect of an enhanced wind of short duration due to a recent TP 
(or pulses). 



can occur if the relative velocity of the AGB wind with respect 
to the ISM is large enough. The thermal emission of the dust in 
the density enhancement at the shock interface between the stel- 
lar wind and the the ISM can be detected in the far-IR (Ueta et 
al. 2006, 2009). While a detached shell produced by TPs would 
be spherical, in the case of wind-ISM interaction the shape of 
the shell will depend on the space motion of the star through the 
ISM. A wind-ISM shell can look spherical as seen for R Cas 
(Ueta et al. 2009) if most of the space motion of the star relative 
to the ISM is in the radial direction. For stars with a high space 
motion (relative to the ISM), the shape of the bow shock will be 
more parabolic with the apex of the parabola in the direction of 
the star's motion relative to the ISM. In the case of CW Leo, the 
shape of the extended emission and the position of the star sug- 
gest that the stellar wind has driven a shock into the ISM. The 
far-IR emission is probably caused by the thermal emission of 
the piled-up dust at the shock interface. 

One may think that the observed UV emission is produced by 
dust scattering of the interstellar light. However, the brightness 
ratio of FUV to NUV is ~10 (see Table[T]i, which is much higher 
than expected in that case (of the order of ~2.4; SC). The only 
other AGB star with UV data probing a wind-ISM interaction is 
Mira (Martin et al. 2007). Martin et al. suggest that collisional 
excitation of cool H2 by hot electrons from the post-shock gas is 
responsible for the UV emission. The faint NUV emission can be 
explained by the H2 emitting in the FUV band. SC suggest that 
the same mechanism may also be the dominant contributor to the 
FUV ring emission in CW Leo. The dust and the FUV emission 
have a similar spatial scale. For the planetary nebula NGC 6720, 
we note that van Hoof et al. (2010) find on the basis of PACS 
and SPIRE data that dust and H2 are co-spatial and argue that 
H2 has been formed on grain surfaces. 

The structure seen in our data may also be the result of both 
mechanisms i.e., mass-loss variation and wind-ISM interaction. 
By estimating the interpulse period for CW Leo, we can check 



Table 1. Derived fluxes for the extended emission 



(yum) 



Flux 

(Jy) 



A 
(pm) 



Flux 

(Jy) 



A 
(jum) 



Flux 

(Jy) 



0.15 
0.23 



(20: 

(2 4 



: 0.2) X 10- 
0.2) X 10-" 



70 
160 



3.51±0.54 
7.70+1.17 



250 4.31±0.66 
350 1.70+0.27 
550 0.73+0.14 



4.2. Space motion of the star and the ISM flow velocity 

From the adopted radial velocity, distance, and proper motion, 
the heliocentric space velocity of the star is derived to be about 
25.1 km s' at a heliocentric inclination angle of the space mo- 
tion vector of 47.8 degrees (measured from the plane of the sky, 
away from us) and a PA of 8 1 .3° for the proper-motion vector in 
the plane of the sky. 

In the PACS/SPIRE wavebands, the observed far-IR surface 
brightness is expected to be proportional to the column density 
of the dusty material in the shell because the optical depth of 
the shell is much lower than unity. Because the bow shock in- 
terface is a parabolic surface arbitrarily oriented in space, the 
column density tends to reach its highest value where the bow 
shock cone intersects with the plane of the sky including the 
central star. Thus, the apparent shape of the bow shock is the 
conic section of the bow shock. Therefore, given the analytical 
formula of the bow shock structure (Wilkin 1996), we can fit the 
observed surface brightness of the bow shock to the conic sec- 
tion of the bow shock cone to derive the heliocentric orientation 
of the bow shock (e.g. Ueta et al. 2008; 2009). From this fitting, 
we determined that the apex of the bow shock cone is oriented 
at 61°9 ± 0°3 (this is degenerate, in a sense that it could be point- 
ing away from us or towards us) with respect to the plane of the 
sky into the PA of 88° with the deprojected stand-off distance of 
(8.0 ± 0.3) X 10'^ cm. 
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Fig. 3. Modified black body (r= 25±3 K and j8= 1.6+0.4) fit 
(dashed curve) to the derived fluxes (symbols) for the extended 
emission. 



Using these numbers in the ram pressure balance equation 
V-i, - ^l -. — — Ta for a mean nucleus number per hydrogen 

Y 4h- //H MH niSM Ri f J b 

nucleus of //h = 1-4, the relative velocity of the star with respect 
to the ISM is V* = (106.6 ± 8.7)/ y/ttisu km s"', where hism is 
the number density of the ISM local to CW Leo in cm"^. 

By following the scheme of Johnson & Soderblom (1987, 
also see Ueta et al. 2008), the heliocentric space velocity com- 
ponents of the star can be converted to the heliocentric Galactic 
space velocity components [U, V, W] of [21.6±3.9, 12.6±3.5, - 
1.8±3.3] km s"' , and also to the LSR Galactic space velocity 
components of [30.6±3.9, 24.6±3.5, 5.2±3.3] km s'K The he- 
liocentric ISM flow velocity is 1 17.6 km s"' if the bow cone is 
facing us (i.e. the apex pointing toward) or 82.6 km s"' if the 
bow cone is facing away from us. 

Wareing et al. successfully modelled bow shocks for AGB 
stars using 3-D hydrodynamical models (RHya and Mira in 
Wareing et al. 2006, 2007a). Their models were able to repro- 
duce in great detail all the components of a bow shock (as- 
tropause, astrotail, vortices) and constrain the space velocity of 
the star, the ISM density, and the mass-loss rate. The models 
were built by varying the velocity of the star relative to the ISM, 
the ISM density and the mass loss rate and the effect of those 
physical parameters on the shape of the bow shock can be seen 
in Wareing et al. (2007b). From the Wareing models, we can 
see that for ISM densities <2 cm"^ (which implies V* >75 km 
s ') and after an evolution of 50 000 yr into the AGB phase, 
the bow shock does not resemble our data. The models have 
a more flattened bow shock with the star closer to the apex of 
the shock. This is for a mass-loss rate between 5x10"^ M© yr"' 
and 10"*" Moyr"', which is an order of magnitude lower than 
the value assumed for CWLeo. For a higher mass-loss rate, we 
would expect an even more important departure from sphericity 
with a significant tail of ejecta. This suggests that the density of 
the local ISM to CW Leo is probably higher than 2 cm"^ imply- 
ing an upper limit of 75 km s^'to V*, the stellar velocity relative 
to the ISM. 
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